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Abstract 

We present a spectral study of the evolution matrix of the 
totally asymmetric exclusion process on a ring at half filling. 
The natural symmetries (translation, charge conjugation com- 
bined with reflection) predict only two fold degeneracies. How- 
ever, we have found that degeneracies of higher order also exist 
and, as the system size increases, higher and higher orders appear. 
These degeneracies become generic in the limit of very large sys- 
tems. This behaviour can be explained by the Bethe Ansatz and 
suggests the presence of hidden symmetries in the model. 

Keywords: ASEP, Markov matrix, symmetries, spectral degen- 
eracies, Bethe Ansatz. 



1 Introduction 

The asymmetric simple exclusion process (ASEP) that plays the role of a 
paradigm in non-equilibrium statistical physics is a model of driven diffusive 
particles on a lattice interacting through hard-core exclusion (Katz, Lebowitz 
and Spohn 1984). Introduced originally as a model for protein synthesis on 
RNA (MacDonald and Gibbs 1969) and for hopping conductivity (Richards 
1977) in solid electrolytes, the ASEP has found a wide spectrum of appli- 
cations ranging from traffic flow (Chowdhury et al. 2000) to surface growth 
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(Halpin-Healy and Zhang 1995, Krug 1997) and sequence alignment (Bund- 
schuh 2002). The ASEP has also motivated many theoretical studies, provid- 
ing a 'minimal' model sophisticated enough to yield a rich phenomenology 
but elementary enough to allow exact solutions [for a review see (Derrida 
1998) and (Schutz 2001)]. 

The ASEP is a memoryless stochastic process governed by simple local 
update rules that can be embodied in a Markov evolution matrix M. This 
matrix contains all the dynamical information about the system: for exam- 
ple, the eigenvector of M with eigenvalue represents the steady state of the 
ASEP and embodies the stationary correlation functions. In particular, the 
matrix product representation of the stationary state (Derrida et al. 1993) 
allows to calculate steady state properties and large deviation functionals 
(Derrida et al. 2003). More generally, the spectrum of M encodes the re- 
laxation behaviour of the ASEP and all non-stationary correlation functions. 
The Markov matrix M can be diagonalized with the Bethe Ansatz that al- 
lows to calculate the energy gap and the dynamical exponent of the ASEP 
(Dhar 1987, Gwa and Spohn 1992, Kim 1995, Golinelli and Mallick 2004). 

Although the steady state of the ASEP and the lowest excitations have 
been extensively studied, little effort has been devoted to investigate global 
spectral properties of the Markov matrix (Bilstein and Wehefritz 1997, Dud- 
zihski and Schutz 2000, Nagy et al. 2002). In this work, we focus on the 
degeneracies in the spectrum of the Markov matrix for the totally asymmetric 
simple exclusion process (TASEP) on a periodic ring at half filling (i.e., the 
density of particles is one-half). The exclusion process has some natural 
invariance properties by translation, charge conjugation and reflection. At 
half filling these symmetries imply the existence of doublets in the spectrum 
of M. However, we find that multiplets with degeneracies of higher order also 
exist for the TASEP. Moreover, the number of these multiplets and their 
degeneracies are given by simple combinatorial factors. These remarkable 
spectral properties suggest that the TASEP at half filling possesses some 
hidden symmetries other than the obvious ones. 

This article is organized as follows. In section 2, we recall the definition 
of the ASEP on a periodic ring and discuss the translation, the charge con- 
jugation and the reflection symmetries, and their effect on the spectrum of 
the Markov matrix M. In section 3, we present some numerical evidence for 
spectral degeneracies in the TASEP at half filling. We show that the com- 
binatorial factors involved can be understood, heuristically, from the Bethe 
Ansatz. In Section 4, we present some concluding remarks and discuss spec- 
tral degeneracies in closely related models. In the Appendix, we analyze the 
natural symmetries of the ASEP with a tagged particle. 
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2 Natural symmetries of the ASEP on a ring 
2.1 The model 

We consider the exclusion process on a periodic one dimensional lattice with 
L sites (sites i and L + i are identical) and n particles. A lattice site cannot 
be occupied by more than one particle {exclusion rule). The state of a site 
i (1 < i < L) can thus be characterized by the Boolean number Tj = 0, 1 
according as the site i is empty or occupied. 

The system evolves with time according to the following stochastic rule: 
a particle on a site i at time t jumps, in the interval between t and t + dt, with 
probability p dt to the neighbouring site % + 1 if this site is empty and with 
probability q dt to the site z — 1 if this site is empty. The jump rates p and 
q are normalized such that p + q — 1. In the totally asymmetric exclusion 
process (TASEP) the jumps are totally biased in one direction (p — 1 and 
'I 0). 

The number n of particles is conserved by the dynamics. The total 
number of configurations for n particles on a ring with L sites is given by 
Q = L\/[n\(L — n)\]. A configuration C can be represented by the sequence 
(ti, t 2 , . . . , t l ). We call ^ t (C) the probability of configuration C at time t. As 
the exclusion process is a continuous-time Markov process, the time evolution 
of iptiC) is determined by the master equation 

^(0 = p(c ) CM(C)- (i) 

The Markov matrix M encodes the dynamics of the exclusion process: the 
element M(C,C) is the transition rate from configuration C to C and the 
diagonal term M(C, C) — — J2c M(C, C) represents the exit rate from con- 
figuration C. 

A right eigenvector vb is associated with the eigenvalue E of M if 

Mijj = Erb. (2) 

The matrix M is a real non-symmetric matrix and, therefore, its eigen- 
values (and eigenvectors) are either real numbers or complex conjugate pairs. 
The spectrum of M contains the eigenvalue E — and the associated right 
eigenvector is the stationary state. For the ASEP on a ring the steady state 
is uniform: the stationary probability of any configuration C is given by 
ip(C) — 1 / O. Because the dynamics is ergodic (i.e., M is irreducible and ape- 
riodic), the Perron- Frobenius theorem (see, for example, Gantmacher 1964) 
implies that is a non-degenerate eigenvalue and that all other eigenvalues 
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E have a strictly negative real part; the relaxation time of the corresponding 
eigenmode is r = —l/Re(E). 

Remark: A configuration can also be characterized by the positions of the 
n particles on the ring, (x\,X2, • • • ,x n ) with 1 < X\ < x<i < . . . < x n < L. 
With this representation, the eigenvalue equation (j2J) becomes 

EiJ)(xi, ...,x n ) = 

^Tp[l/}(xi,...,Xi-i, Xi-1, x i+ x, . . . ,x n ) -tp{xi, . . . ,x n )] + 
i 

q bPi^i, • • • , Xj-x, xj + 1, x j+1 , ...,x n )- ip(xx, x n )) , (3) 

j 

where the sum runs over the indexes i such that < Xj — 1 and over the 
indexes j such that x } ■ + 1 < Xj + \ , i.e., such that the corresponding jumps 
are allowed. 

2.2 Translation invariance 

The translation operator T shifts simultaneously all the particles one site 
forward: 

T(ri, r 2 , . . . , t l ) = (r L , n, r 2 , . . . , tl-i) ■ (4) 

Because of periodicity we have T L = 1. Thus, the eigenvalues (or impulsions) 
k of T are simply the L-th roots of unity: 

k = e 2iTm/L for m = 0, . . . , L - 1 . (5) 

We denote by T(k) the eigenspace of T corresponding to the impulsion k. 
The projection operator Tk over this eigenspace is given by 

r fc = )p-r. (6) 

L r=0 

The complex conjugation transforms a T-eigenvector with eigenvalue k into 
an eigenvector with eigenvalue k* = k~ l , i.e., T(k*) = [T(fc)]*. 
The ASEP on a periodic lattice is translation invariant, i.e., 

MT = TM . (7) 

The matrix M and the translation operator T can therefore be simultaneously 
diagonalized: M leaves each eigenspace T(k) invariant. We denote by Sp(fe) 
the set of the eigenvalues of M restricted to T(k). Using complex conjugation 
we obtain the property 

Sp(F) = [Sp(*)r • (8) 
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2.3 Right-Left Reflection 

The reflection operator R reverses the right and the left and is defined by 

R(n,T 2 , ...,t l ) = (tl,t l _i, . . . ,T 2 ,Ti) . (9) 

We have R 2 — 1; the eigenvalues of R are thus r = ±1. The reflection i? 
reverses the translations, i.e., 

RT = T~ 1 R, (10) 

and transforms a T-eigenvector with eigenvalue k into an eigenvector with 
eigenvalue k* = k. The operator R does not commute with the Markov 
matrix M because the asymmetric jump rates (p ^ q) are not invariant 
under the exchange of right and left. Writing explicitly the dependence of 
M on the jump rates, we have 

RM^R- 1 = M (q , p) . (11) 

Thus, in general the reflection is not a symmetry of the ASEP (it is only a 
symmetry of the symmetric exclusion process). 

2.4 Charge conjugation 

The charge conjugation operator C exchanges particles and holes in the sys- 
tem, i.e., a configuration with n particles is mapped into a configuration with 
L — n particles: 

C(n, T 2 , . . . , T L ) = (1 - Ti, 1 - T 2 , . . . , 1 - T L ) . (12) 

The operator C satisfies the relations: 

C 2 = l, CT = TC, CR = RC. (13) 

By charge conjugation, particles jumping forward are mapped into holes 
jumping forward. But holes jumping forward are equivalent to particles 
jumping backward. Writing explicitly the dependence of M on the number 
of particles and on the jump rates, we have 

CM^C- 1 = M {L _ niq , p) . (14) 

We notice that the number of particles is conserved by C only at half filling 
[L = 2n). Thus, the charge conjugation is not a symmetry of the ASEP 
except for the symmetric (p = q) exclusion process at half filling. 
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2.5 CR symmetry 

Using equations fTTf and (JHJ) to combine the charge conjugation C with the 
reflection R, we obtain 

{CR)M {n ^ q) {CRY l = M (L _ n , Pi5) . (15) 

Thus, for a given L, the CR operator maps the ASEP with n particles into the 
ASEP with the same jumping rates but with L — n particles. This implies 
that the spectrum of M for n particles is identical with the spectrum for 
L — n particles (because CR transforms eigenvectors of M n into eigenvectors 
ofM L _ n ). 

Hereafter, we shall always consider the ASEP model at half filling, i.e., 
L = In. In this case, the CR operator constitutes an exact symmetry because 
(CR)M = M(CR). 

The ASEP at half filling is therefore invariant under each of the two 
symmetries, translation T and CR. Note that they do not commute with 
each other. Rather, we obtain from equations (jlOjl and (|13jl 

(CR)T = T-\CR). (16) 

Hence, the CR transformation maps the subspace T(k) into T(k*) and thus 

Sp(*) = Sp(F) . (17) 

For k 7^ ±1, the two subspaces T(A;) and T(fc*) are distinct and equation (fI7|l 
implies that the corresponding eigenvalues of M are doubly degenerate. For 
k = ±1, the CR transformation leaves T(k) invariant. Thus, the natural 
symmetries do not predict any degeneracies in the sets Sp(±l). However, 
each of the subspaces T(±l) is split into two smaller subspaces invariant 
under CR and on which CR = ±1. 

Equations (jBJ) and (fTTj) imply that for all k the set Sp(fc) is self-conjugate, 

= [s P (*)r . (is) 

This means that Sp(fc), for all k, is made only of real numbers or complex 
conjugate pairs. 

In the Appendix, we explain how a charge conjugation operator C can be 
defined for the ASEP with a tagged particle so that CR remains a symmetry. 

3 Spectral degeneracies in the TASEP 

The above discussion, based on obvious symmetries of the ASEP at half 
filling with generic jump rates p and q, suggests that the spectrum of the 
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Markov matrix should be composed of singlets (for impulsion k = ±1) and 
doublets (for k ^ ±1); so multiplets of higher order have no reason to appear 
and therefore should not exist generically. 

In this section we describe spectra for systems of size L < 18 at half filling 
obtained from numerical diagonalization for the TASEP, i.e. the model where 
jumps are allowed only in one direction (q = 0). These spectra contain many 
unexpected degeneracies of very special orders. Our numerical observations 
lead us to conjecture some formulae for general L. We shall show that these 
formulae can be derived, at least heuristically, from the Bethe Ansatz equa- 
tions. The existence of such degeneracies of higher orders strongly suggests 
the presence of unknown underlying symmetries. 

3.1 Numerical results and conjectures 

The spectral degeneracies of the Markov matrix M are numerically calculated 
with the following procedure. First, the matrix M is built for a given size L 
and number of particles n. As the system is periodic, we use the translation 
invariance T in order to split the space of configurations into L subspaces 
T(k) with impulsion k (see section Thus M is block-diagonal with L 

blocks. The spectrum of each block is then computed by a routine of the 
Lapack library (Anderson et al., 1999) with the standard numerical double 
precision (around 15 decimal digits). The L spectra obtained by numerical 
diagonalization are then collected together in order to search and count the 
degeneracies. 

The time required to compute the complete spectrum of a matrix of size 
iV is proportional to N 3 . Since the sizes of the blocks T(k) grow as 2 L , the 
computation time grows as 8 L . Thus only small systems can be studied by 
this method: we studied systems of size L < 18. 

The problem is now to choose a criterion to decide when two eigenvalues, 
computed with finite precision, are equal. By calculating the differences 
between all the eigenvalues we observe that these differences are always either 
smaller than 10 -12 or larger than 10~ 4 . Consequently we deduce that two 
eigenvalues with a difference smaller than 10 -12 are in fact equal and that this 
apparent discrepancy is due to the finite precision of the computation. This 
allows us to enumerate the degeneracies. We emphasize that the existence of 
this gap of eight orders of magnitude (between 10~ 12 and 10~ 4 ) in the data 
of the differences is crucial. If the distribution of the differences was more 
continuous (i.e., if there was no obvious gap), the choice of a criterion to 
decide which eigenvalues are identical would have been much more difficult 
and somehow arbitrary. 

Let m(d) be the number of multiplets with degeneracy of order d, i.e., the 
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L 



d=l d = 2 d = Q d = 20 d = 70 



2 
4 
6 

8 

10 
12 
14 
16 
18 



3432 
12870 
48620 



2 
6 
20 
70 
252 
924 



2 

4 1 

8 6 

16 24 1 

32 80 10 

64 240 60 1 

128 672 280 14 

256 1792 1120 112 1 

512 4608 4032 672 18 



Table 1: Spectral degeneracies in the TASEP at half filling: L is the size of 
the lattice; Q is the dimension of the configuration space; the other columns 
give m(d) the number of multiplets of degeneracy d. These numbers are given 
by the simple combinatorial formulae 17^) and (fl/^). 

number of sets of d equal eigenvalues. By definition, we have J2d d m{d) = £1, 
the total number of eigenvalues. For the TASEP model with half-filling, the 
numbers m(d) are given in Tabled for sizes L < 18. Although only singlets 
(d = 1) and doublets (d = 2) are predicted by the natural symmetries T 
and CR, numerous degeneracies of higher order also appear. We notice from 
Table [U that the observed orders of degeneracies are only 1,2,6,20,70. These 
numbers obey the sequence 



for a system of size 2n, r takes all integral values from to n/2. From an 
empirical study of Tabled we also remark that for L = 2n < 18, 



For r = and 1, these equations provide the number of singlets m(l) = 2 n 
and the number of doublets m(2) = n(n — l)2 n ~ 3 . We conjecture that equa- 
tions ()19|) and (|20|) are true for all L. We have verified that these relations 
satisfy the sum rule J2d d m(d) = Q: 




(19) 




(20) 



r>0 0<r<n/2 



d r m(d r ) = 



(n — 2r)\ r! r! 



2n 



n 



) 



(21) 



The last equality is obtained by identifying the terms of order x° in the 
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identity 

v x/ ni+n. 2 +n 3 =n n 1 [n 2 ln 3 l X k \ / 

Using equations (fT9"j) and (J2U|l . we can analyze the large L (or n) be- 
haviour. For example, we note that the number of singlets, 2™, becomes 
negligible with respect to the total number of eigenvalues, f2 ~ 4 n /n 1//2 , as 
L grows. As explained in section 12. 5[ the CR symmetry implies that the 
singlets necessarily belong to the subspaces T(±l). However, the number 
of these singlets is much smaller than the dimension of subspaces T(±l) 
(which is approximately given by 0,/n oc A n /n 3 ^ 2 ). Therefore, multiplets 
must necessarily exist in the subspaces T(±l): this fact is an indication for 
a hidden symmetry that would commute with T but not commute with CR, 
for example. 

More precisely, the number m(d) of multiplets is maximal for r ~ n/6 
which corresponds to a degeneracy of order d oc 2 ra / 3 with m(d) oc 3 n . Then, 
we obtain that dm(d) oc (2 1 / 3 3) n < Q because 2 1 / 3 3 « 3.78 < 4. On the 
other hand, the number dm(d) of eigenvalues counted with their degeneracies 
is maximal for r ~ n/4. In this case, we have d oc 2™/ 2 , m(d) oc 8 n / 2 and 
dm(d) oc 4 n oc Q. Loosely speaking, we deduce from this estimation that 
most of the eigenvalues in the spectrum are degenerate with a degeneracy of 
order T 1 ! 2 . 

3.2 Bethe Ansatz analysis of the spectral degeneracies 

The Markov matrix of the exclusion process can be diagonalized thanks to 
the Bethe Ansatz (Dhar, 1987) that assumes that an eigenvector ip of M can 
be written as 

i/j(xi, . . . , x n ) = ^2 A a z^ z*( 2 ) . . . z%fa , (23) 

where S n is the group of the n\ permutations of n indexes. The coefficients 
{A a } and the wave-numbers {zi, . . . , z n } are complex numbers to be deter- 
mined. For TASEP at half filling, it is convenient to use the fugacity variables 
Zi = 2/zi — 1 and to introduce an auxiliary complex variable Y. The Bethe 
equations can then be reformulated as explained below. (For more details, 
see Gwa and Spohn 1992; Golinelli and Mallick 2004). In terms of the n-th 
roots of Y defined as 

y m = yi/n e (m-i)2Wn for m = 1, . . . , n , (24) 



Golinelli, Mallick — Hidden symmetries in ASEP 



10 



we calculate 2n numbers (Zi, . . . , Z 2n ) given by 

Z m — (1 — Vm) 1 ^ 2 ] Z m+n = —Z m . (25) 

In order to select n fugacities among (Zi, . . . , Z 2v ), we introduce a choice 
function c : {1, . . . , n} — » {1,..., 2n} that satisfies 1 < c(l) < . . . < c(n) < 
2n . The Bethe equations now become equivalent to the self-consistency equa- 
tion 

MX) = Y, (26) 
with A C {Y) = (27) 

j=l Z c{j) + 1 

Given the choice function c and a solution Y of this equation, the Z c ^s are 
determined from Eq. (|2*5j) and the corresponding eigenvalue E c is given by 

n 

2E C = -n + J2 ZcU) ■ (28) 

3=1 

We now assume that the Bethe Ansatz is complete, i.e., the Bethe equa- 
tions yield a complete basis of eigenvectors for the ASEP. Besides, we also 
note that the number of different choice functions is Q = (2n)!/n! 2 , Q be- 
ing precisely the size of the Markov matrix. We are thus led to make the 
stronger hypothesis that for each choice function c (among the Q possible 
choice functions), the self-consistency equation (J26|) has a unique solution Y 
that yields one eigenvector ip c and one eigenvalue E c . We further assume 
that these eigenvectors are linearly independent. 

Using this one-to-one hypothesis, we will now explain the spectral degen- 
eracies and derive equations (fT^j) and (|2~0"j). We recall that the set (Zi, . . . , Z 2n ) 
is constituted of n pairs, each pair containing two opposite fugacities (Z m+n = 
—Z m ). A choice function c selects n numbers in the set (Zi, . . . , Z 2n ); amongst 
these n selected fugacities, 2r appear in pairs and n — 2r appear as singles, 
with < 2r < n. The crucial remark is that for a given choice function c the 
Bethe equations and the eigenvalue E c depend only on the n — 2r selected 
single fugacities. This fact is clear from equations (J2TH 1271 and I2EJ) because 
contributions of Z m and Z m+n cancel each other. Thus two different choice 
functions c\ and c 2 that contain the same single fugacities (but different 
pairs) lead to the same eigenvalue but to different eigenvectors. Hence all 
the choice functions that contain the same single fugacities form a multiplet 
with the same eigenvalue. The number of these choice functions is, thus, the 
order of degeneracy of the corresponding multiplet. This order of degeneracy 
is determined solely by the integer r. 



Golinelli, Mallick — Hidden symmetries in ASEP 



11 



For a given value of r (with < 2r < n), we now calculate (i) the 
number d r of choice functions in each multiplet, i.e., its degeneracy and (ii) 
the number m(d r ) of different multiplets. 

(i) Two choice functions c and d are in the same multiplet if they differ 
only in the choice of the r pairs. The n—2r single fugacities having been fixed, 
we must choose r pairs amongst 2r available pairs. The number of choice 

functions in a multiplet is therefore given by d r = ^ ^ ^ in agreement with 

equation (JTHJ). 

(ii) A multiplet is characterized by its s = n — 2r single fugacities. The 
number of different ways to select these single fugacities is given by 



in agreement with equation ((201) • 

These calculations thus explain the numerical results described in the 
previous section. However, our analysis is based on the hypothesis of a one- 
to-one correspondence between choice functions and linearly independent 
eigenvectors. We emphasize that this one-to-one hypothesis is more strin- 
gent than the assumption of the completeness of the Bethe Ansatz. At the 
present stage, we do not have a proof of this hypothesis but we have tested 
it numerically for small system sizes: it is verified for L — 2,4, 6, 8. 

4 Conclusion and Generalizations 

In this work, we have studied the degeneracies in the spectrum of the Markov 
matrix of the totally asymmetric exclusion process on a periodic ring at half 
filling. Numerical results show the presence of many degeneracies with spe- 
cific orders. These orders and the number of multiplets of a given degeneracy 
can be calculated from the Bethe Ansatz equations under the assumption of 
a one-to-one correspondence between the eigenvectors of the Markov matrix 
and some choice functions of the fugacities that appear in the Bethe Ansatz. 
Certainly, a proof (or a disproof) of this hypothesis would be of interest. We 
believe that this assumption must be true, at least in a weaker form, because 
it leads to a correct enumeration of the spectral degeneracies and multiplets. 

We emphasize that the observed higher order degeneracies are not pre- 
dicted by the natural symmetries of the TASEP. Indeed, we showed that 
invariance under translation and charge conjugation plus reflection, predicts 
only the existence of two fold degeneracies in the spectrum. Therefore, we 



m(d r ) 



2n{2n - 2) . . . (2n - 2{s - 1)) 
7\ 
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believe that the ASEP presents some hidden symmetries, that are responsi- 
ble for the numerous and large generic degeneracies in the spectrum of the 
Markov matrix. By inspecting numerically the multiplets, we have observed 
that each of them are composed of eigenvalues with different impulsion k: 
thus a hidden symmetry does not necessarily commute with the translation 
operator T. A better understanding of these symmetries would be of great 
significance. 

We have also analyzed the degeneracies in some related models in order 
to study the influence of asymmetry and periodicity. For the partially asym- 
metric exclusion process with non zero left and right jump rates (p > q > 0) 
on a periodic lattice and at half filling, numerical computations show the 
existence of singlets and doublets and the absence of degeneracies of higher 
order. However, the number of singlets is exactly 2™ (i.e., the same as for 
the TASEP). From the natural symmetries, one would expect a much larger 
number of singlets (of the order of 4 n /n 3 / 2 ). This is again a hint for some 
hidden symmetries in the partially asymmetric exclusion process (Alcaraz et 
al. 1994). The TASEP can be also defined on an open lattice: it presents 
a rich phenomenological behaviour that was thoroughly investigated using 
the matrix product representation (Derrida et al. 1993). In this system the 
translation symmetry T is lost and the number of particles is not conserved. 
We have numerically computed spectra for different sizes L and different en- 
trance and exit rates: no degeneracies have been observed in any of these 
cases. 

Another direction of interest is to study the TASEP model for arbitrary 
filling. In fact, we have numerically observed degeneracies of high order in 
the spectra when the density is a simple fraction. Degeneracies at arbitrary 
filling can be analyzed by Bethe Ansatz but the calculations seem to be 
far more complicated. General results will be presented in a future work 
(Golinelli and Mallick 2004a). 
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Appendix : Natural symmetries of the ASEP 
with a tagged particle 

To study the motion of particles in the ASEP, it is often interesting to tag 
one of them (Spohn, 1991). All particles are then automatically labelled 
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because overtaking is forbidden. A configuration C of n particles on a ring 
with L sites can now be specified by the sequence (ti, t 2 , . . . , Tl', m) where 
Tj = 0, 1 according as the site i is empty or occupied and the tag is on 
the m-th occupied site with 1 < m < n (m is defined modulo n). The 
number of configurations of the ASEP with a tagged particle is given by 
nVL = nL\/[n\(L-n)\}. 

The translation operator T shifts simultaneously all the particles one site 
forward, i.e. 

T(n, r 2 , . . . , r L ; m) = (r L , n, . . . , r L -\\ m + r L ). (30) 

Indeed, the index m increases by 1 if and only if T moves a particle from site 
L to site 1. The ASEP with a tagged particle remains translation invariant: 
let M be its Markov matrix, we have MT = TM. As T L = 1, the eigenvalues 
k of T are the L-th roots of unity [see Eq. (J5J)]. 

The ASEP with a tagged particle presents a new symmetry P, which 
consists in shifting the tag one particle forward but without moving the 
particles, i.e. 

P(ri, . . . , r L ; m) = (n, . . . , r L ; m + 1) . (31) 

As P n = 1, the eigenvalues a of the tag-shift operator P are simply the n-th 
roots of unity. We have 

PM = MP; TP = PT . (32) 

The matrix M can thus be simultaneously diagonalized with T and P. We 
call E(A;, a) the common eigenspace of T and P corresponding to the eigen- 
values k and a respectively, where 

k = e 2inj/L for j = 0, . . . , L - 1 , (33) 
a = e 2ina/n for a = 0, . . . ,n - 1 . (34) 

and Sp(/c, a) the set of the eigenvalues of M restricted to E(fc, a). 

We remark that the Markov matrix of the ASEP without a tag is the 
restriction of the Markov matrix of the ASEP with a tag to the subspace 
a — 1. Indeed, the projection on the subspace a = 1 of a configuration 
(ti, r 2 , . . . , tl; m) is given by S" =1 (ri, r 2 , . . . , r L ;j)/n; this projection renders 
the particles indistinguishable. 

Let us now consider a configuration (T 1; r 2 , . . . ,t l ; m), z.e., the tagged 
particle is the m-th particle from the left. With the left-right reflection R, 
it becomes the m-th one from the right, so the (n + 1 — m)-th one from the 
left, i.e., 

R(n, r 2 , . . . , r L ; m) = (r L , . . . , r 2 , n; n+l-m) . (35) 
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We recall that R is not a symmetry of the ASEP because the right and left 
jump rates, p and q, are swapped. 

The action of the charge conjugation operator C is more subtle. The op- 
erator C transforms a system with n particles (including one tagged particle) 
into one with L — n particles (including one tagged particle). The dimen- 
sions of the corresponding configuration spaces will be the same only in the 
half filling case; thus for the tagged ASEP, the charge conjugation C can be 
well defined only for L = In. Moreover, because C transforms the tagged 
particle into a hole, the tag cannot remain on the same site and we must 
specify where it goes. One solution, that will ensure the CR-invariance of 
the ASEP, is to define C as follows: 

C(n, r 2 , . . . , t l \ m) = (1 - n, 1 - r 2 , . . . , 1 - r L ] n+l-m). (36) 

Hence, after charge conjugation, the tag is on the m-th particle from the 
right. Consequently, we obtain, at half filling 

CR(n, t 2 , . . . , t l ; m) = (1 - r L , ■ ■ ■ , 1 - r 2 , 1 - n; m). (37) 

It follows that 

(CR) 2 = 1 ; P(CR) = (CR)P . (38) 
It is also possible to show that CR is a symmetry of the model: 

M(CR) = (CR)M . (39) 

(a little care is needed when analyzing jumps between sites L and 1). 

Combining Eq. (J3*U|) with Eq. ([37)1 . the relation between CR and T is 
obtained 

CRT(t 1 ,t 2 ,. . . ,r L ;m) = (1 - Tfc_i, . . . , 1 - r x , 1 - r L \ m + t l ) , (40) 
TCRT(t 1 ,t 2 ,. . . ,r L ;m) = (1 - r L , 1 - t l _ x , . . . , 1 - n; m + 1) , (41) 

so TCRT = PCR, or equivalently 

T (CR) = (CR) PT- 1 . (42) 

Hence the CR transformation maps the subspace E(fc, a) into E(ak*, a) and 
thus 

Sp(k, a) = Sp(ak*,a). (43) 

For k 2 ^ a, these two subspaces are distinct but have the same spectrum: 
their eigenvalues are doubly degenerate. For k 2 = a, the CR symmetry 
leaves E(fc, a) invariant and splits it into two smaller subspaces with cr = ±1, 
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without degeneracy. With the notations of Eqs. (jHHllMj) . the condition k 2 = a 
is equivalent to j = a or a + n. 

Complex conjugation implies that [Sp(fc, a)}* = Sp(k*,a*). Therefore, 
we obtain from Eq. (|4*3*|) that Sp(fc, a) is self-conjugate if a = 1 for any k 
(this case is equivalent to the ASEP without tag) or if a = — 1 for k = ±1. 

To summarize, the CR symmetry can still be defined at half filling for 
the generic ASEP with a tagged particle. This symmetry together with 
translation invariance predicts the existence of doublets in the spectrum [see 
Eq. (@3D]. 

For the totally asymmetric model (TASEP) with a tagged particle, nu- 
merical diagonalization shows numerous degeneracies of order higher than 2 
inside each P-invariant subspace. It follows that hidden symmetries should 
also exist in this system and should commute with the tag-shift operator P. 
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